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Abstract — Spatially-coupled (SC) codes are constructed by 
coupling many regular low-density parity-check codes in a chain. 
The decoding chain of SC codes stops when facing burst erasures. 
This problem can not be overcome by increasing coupling 
number. In this paper, we introduce multi-dimensional (MD) 
SC codes. Numerical results show that 2D-SC codes are more 
robust to the burst erasures than ID-SC codes. Furthermore, we 
consider designing MD-SC codes with smaller rateloss. 

I. Introduction 

Spatially-coupled (SC) low-density parity-check (LDPC) 
codes attract much attention due to their capacity-achieving 
performance and a memory-efficient sliding-window decoding 
algorithm. The studies on SC-LDPC codes date back to the 
invention of convolutional LDPC codes by Felstrom and 
Zigangirov [T]. Lentmaier et al. observed that (4,8)-regular 
convolutional LDPC codes exhibited the decoding perfor- 
mance surpassing the belief propagation (BP) threshold of 
(4,8)-regular block LDPC codes \2\. Further, the BP threshold 
coincides with the maximum a posterior (MAP) threshold of 
the underlying block LDPC codes with a lot of accuracy. 
Constructing convolutional LDPC codes from a block LDPC 
code improves the BP threshold up to the MAP threshold of 
the underlying codes. 

Kudekar et al. named this phenomenon "threshold sat- 
uration" and proved rigorously for the binary-input erasure 
channel (BEC) [31 and the binary-input memoryless output- 
symmetric (BMS) channels. lH. In the limit of large di,dr,L 
and w, the SC-LDPC code ensemble {di, dr, Ljw) 111 was 
shown to universally achieve the Shannon limit of the binary- 
input memoryless symmetric-output (BMS) channels under BP 
decoding. 

In this paper, we deal with a serious problem of SC- 
LDPC codes. SC-LDPC codes are constructed by couphng 
L regular LDPC codes in a chain. Belief propagation (BP) 
is employed to decode the chain of codes starting from the 
end points of the chain. The BP decoding of SC codes stops 
when facing burst erasures. In other words, the decoding error 
probability remains positive from the section at which the 
burst erasures are received. This problem can not be solved 
by increasing coupling number In this paper, we introduce 
multi-dimensional (MD) SC codes to overcome this problem. 
Numerical results show that 2D-SC codes are more robust to 
the burst erasures than ID-SC codes. Furthermore, we consider 



designing MD-SC codes with small rateloss as 0(1/ L^), 
where D is the dimension. 

II. Multi-Dimensional Coupled Codes 

A. Definition: {di,dr, L,lu, Z) codes 

Definition 1: Define Zl := Z/LZ = {0, 1, . . . , L - 1}. For 
a,b e Z^. Consider sections on D-dimensional discrete 
torus Z|^. For bit node degree rf/ > 3 and check node degree 
dr > di, a coupling number L > w, connecting rate < cjj < 
1 {j € Z|^), and a shortened domain Z C Z£, we define ML- 
SC {di,dr,L,u!, Z) codes as follows. Throughout this paper, 
we fix utj as 

_ jl/w^ j e[0,w- 1]^ 
1 otherwise, 

where we denoted [a, b] := {a, a + 1, . . . , 6 — 1, 6}. 

Each section i G [0, L — 1]^ has M bit nodes of degree di 
and j[-M check nodes of degree dr. Connect edges between 
bit nodes and check nodes uniformly at random so that bit 
nodes in section i are connected to check nodes in section 
i + j (.? £ ^l) ^i'-h ujjM edges, respectively. Shorten the 
bit nodes in section i&Zc Z|J'. Namely, the shortened bit 
nodes are set to and are not transmitted through the channel. 

Discussion 1: In 111, spatially-coupled codes of coupling 
number L were defined over section [—L, +L] and the bit 
nodes outside [—L, +L] were shortened. Some might think it 
is more natural to define MD-SC codes over [— L, than 
over [0, L — l]^ and shorten the bit nodes outside [— L, 
If we defined so, it would be difficult to distinguish the 
effect of MD extension from the boundary effect from each 
dimension as ID-SC codes. This is why we employ the codes 
in Definition [T] 

Lemma 1: The coding rate 
{di,dr, L, u), Z) codes is given by 



R{di,dr, L,Lu, Z) of 



dr LD - #Z 



(1) 



Proof. We will count the number of transmitted bit nodes 
and valid check nodes. Let V and C denote these numbers, 
respectively. Since check nodes adjacent only to shortened bit 
nodes are not giving any constraint on the code, it is sufficient 



to count the check nodes adjacent to unshortened bit nodes. 
Since the degree of check nodes are dr, a check node in 
section i has dr edges connecting to shortened bit nodes with 
probability 



( E 



Therefore, the average number of check nodes which are 
adjacent to at least one unshortened bit nodes is given by 



C 



di 



j:i+jeZ 



There are V — M{L^ — #Z) unshortened bit nodes. We 
calculate the coding rate as 1 — C /V, which concludes the 
lemma. □ 

III. Density Evolution Analysis 



We consider the transmission takes place over the BEC(e) 

with erasure probability e. The BP decoding is employed. Let 

(I) 

PI denote the erasure probability of BP messages from bit 

nodes to check nodes at the £-th iteration round. Let denote 
the erasure probability of BP messages from check nodes to 
bit nodes at the £-th iteration round. Since the bit nodes in 
section in Z are shortened, pf'' are given as 



(2) 



For t > 1, PI = Q for shortened section i^Z and 



(I) 1 (. (f-i) 



for i^Z. The decoding erasure probability Pj^^^ is given by 



We define the BP threshold e* as 

e* = sup{e > I lim = 0}. 

Namely, for the erasure probability below the threshold e*, the 
decoding erasure probability goes to zero. 

A. Shortening hyperplane of width w 

Choose a hyperplane of width w in D-dimensional space 
is chosen as shortened domain Z. The following proposition 
asserts that the density evolution is equivalent to that of of ID 
system. 

Proposition 1: Let us define 

Z:= [0,u;~l], 

Z := {i= eZf I in e [0,^-1]}. 



We use ~ for variables for ID system throughout this paper for 
the sake of readability. Then we have 

e*{dr, dr,L,Lj,Z) = e*{dr, dr, L,uj, Z), (3) 

R{dr, dr,L, Z) = R{dr,dr, L, LO, Z) (4) 

= l-^-0{l/L). (5) 



Proof: We give a proof for D = 2. The proof for D > 2 
allows similarly. It is si 
> 0. It is obvious that 



follows similarly. It is sufficient to show pl^ — pi for any 



~(o) _ (0) _ 1 41 e [0, w - 1] 

-P(iui2) " |g otherwise, 

for 11,12 e Zl. Assume pf^ — p|^| for i. From Eq. (|2|i 
and the definition of uj and lo, it follows that 

-. (i; — luj — 1 ji 



d,-l 



qf.'''\ = i-[i-A,yypf , 

^(n,*2) \ yj2 Z-^ /-^ ^(il-3l,i-2-32) 

ii=Oi2=o 

^ w — lw — 1 



di-l 



yj^ X ( X ( yi-r J 1,12+32) 

31=0 j2=0 

31=0 

Thus we have n,-^' = p''/:^ . , for any £ > 0, which concludes 
(13]l. We derive ^ as follows. 



R[dr, dj., L, w, Z 

(il,i2)eZ| Ui,32)-{il+3l,i2+32)eZ 



di/dr 
L2 - wL 

di/dr 
L — w 



E E 

^ E E ^0ij2))' 

(i 1 ,0) ezi (ii ,j2 ) : (ii +ii , j2 ) e 2: 

-u.) — 1 

E (i-(E EiM^'- 



(ii.o)ez| i2=o ji:0<ii+ji<u)-i 

= R{dr, dr., L, UJ, Z). 

Equation ^ follows from 

dzi-«^-2Er:o(^)''- 



R{dr,dr, L^Lo, Z) = (1 



L — w 



of which proof is appeared in f3] for coupled codes defined 
on [-L,L\. □ 
In the next section, we will see these MD-SC codes with 
shortening domain as a hyperplain behave differently from the 
ID-SC codes. 
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Fig. 1 . Transition of message en'or probability at each section i and at 
iteration £ of ID-SC (d; = 3, dr = 6, L = 101, 6j,Z = {0, ±1}) codes 
with to = 4. The channel is BEC(e = 0.48) with 2 burst section erasures 
injected. Decoding stops around at section i = ±31 where burst erasures are 
injected as p''^' 
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Fig. 2. The BP threshold of ID-SC and 2D-SC (d; = 3, dr = 6, L, u), Z) 
with Z as hyperplane of width w. The BP threshold of ID-SC codes are 
badly degraded when burst section erasures exits. The 2D-SC codes are more 
robust than the ID-SC codes. 



IV. Robustness for Burst Erasure 

In this section, we consider burst erasures and demonstrate 
robustness of 2D coupled codes. Spatially coupled codes are 
constructed by coupling L regular LDPC codes of length M. 
Assume that we are transmitting bits coded by ID-SC codes 
and a burst section erasure of length M occur at section i. We 
call such burst erasures for a section a burst section erasure. 
Decoding proceeds from the section in Z. Such burst section 
erasure is described as e^^ = = 1. 

Can ID-SC codes correct such burst erasures? Figure [l] 
shows the transition of message error probability of ID-SC 
(rf/ = 3, dr = 6, L = 101, d), Z = {0, ±1}) codes with w = A. 
The channel is BEC(e = 0.48) with 2 burst section erasures 
injected. Decoding stops around at section i ~ ±31 where 
burst erasures are injected as pf^ = 0.52. The ID-SC codes 
can not recover such burst erasures. 

Figure ID shows the transition of decoding error rate of 2D- 
SC {di ^ 3,dr ^ 6,L = 101, w, Z) codes with u)^ = 4 and Z 
as square segment of size 15. The channel is BEC(e = 0.48) 
with 20 burst section erasures injected. Each burst section 
erasures are described as pf^ = 1.0. The 2D-SC codes are 
capable of recovering such burst section erasures. 

Figure 12] compares the BP threshold values of ID-SC codes 
and 2D-SC codes with width-u> hyperplane shortening domain. 
The degrees di and dr are set to 3 and 6, respectively. We in- 
jected one or two burst section erasures. The coupling number 
L for each plotted point is chosen sufficiently large so that 
the BP threshold is not increased due to the rateloss, the burst 
error sections are not affected each other, and each plotted BP 
threshold value converges. Note that the BP threshold value 
is about 0.4882 when there is no burst section erasures. For 
small coupling window size w > 4, the BP threshold of ID- 
SC codes is 0. This is badly degraded from 0.4882. This 
degradation would not be mitigated by increasing L. When 



w — 2 the BP threshold is 0, namely the burst section erasure 
is not recovered even if all other sections were recovered. This 
can be explained by the theorem in the next section. 

On the other hand, 2D-SC codes are not degraded from the 
case of no burst section erasures even if w is small. From this 
observation, 2D-SC codes are more robust to burst section 
erasures than ID-SC codes. The x-axis indicates . This 
is intended to be dealt fairly with respect to the number of 
coupled neighboring sections both at ID and 2D. 

A. Bound on Performance 

In the previous section we observed that ID-SC codes did 
not recover one single burst section erasures when w ~ 2. 
This can be explained by the following theorem. 

Theorem 1: The MD {di, dr, L, lu, Z) code of dimension D 
can not recover one single burst section erasure at section i 
if e^_ > e^^{di,dr)w^ if i ^ Z, where e'^^{di,dr) is the BP 
threshold of uncoupled {di,dr) codes. 

Proof: Let us consider the best case, namely other sec- 
tions have no erasures. To be precise, ej = for j ^ i. The 
density evolution equations can be written as 



P 



(j ^ i) 



-1) 
+j-k) 



1 



1 



(f-l)xd,--l 



d,-l 



Denoting p. 



it) 



P. 



PI 



'^^'^ /w^, we have 



This can be viewed as density evolution of uncoupled {di,dr) 



code over BECiei/w^). Hence 
which concludes the theorem. 



> Oif ■ 



>e^^{di,dr). 



From Theorem[T] one can see that {di = 3, dr = 6, L, uj, Z) 
with w = 2 can not recover one single burst section erasure 
since e^^l> 0.4294 x 2 0.8588 = e^^ {di , dr)w^ . The 
BP threshold gets degraded even worse when the number of 
burst section erasures is 2. 

V. Rateloss Problem of Multi-Dimensional SC 
Codes and Its Mitigation 

As one can see in ([T]i, the rate of SC codes is less than 
the uncoupled codes 1 — The ID-SC codes have rateloss 
0{1/L). The ID-SC codes could have rateloss 0(1/ L^) by 
coupling sections as ID-SC codes. The _D-dimensional SC 
codes with hyperplane shortened domain have only 0{1/L) 
while there are sections. This is a problem. Is it possible 
to design MD-SC codes with rateloss 0{1/L^) by keeping 
the BP threshold the same? 

Define the shortening domain as a hypercube of size z as 
follow. 

We claim that the rateloss of the codes with this Z has rateloss 
0(1/ L^). The number C of check nodes that are adjacent to 
unshortened bit nodes is not greater than the number of all 
check nodes. 

C < ^ML"" 

dr 

There are V = M{L^ — z^) unshortened bit nodes. Thus we 
have the coding rate as 



R=l- C/V > 1 



L 



D 



1 



di_ 

dr 



O 



dr - Z^ 



Note that we are not saying that this rate is better than the 
coding rate of ID-SC codes. It is fair to compare the coding 
rate keeping the number of sections Lc the same. From this 
point of view, the rateloss of both ID-SC codes and the MD- 
SC codes scales with 0(1/ Lc), where Lc = L for ID-SC 
codes and Lc = L^ for MD-SC codes of dimension D. 

Does the BP threshold attain the MAP threshold of the 
uncoupled codes? Figure |3] shows the BP threshold of 2D-SC 
{di = 3,dr = 6, L,w = 2) codes with shortening domain Z as 
hypercube of size z. We take sufficiently large coupling num- 
ber L so that each plotted point converges. We observe that 
the BP threshold approach the MAP threshold of uncoupled 
codes as z gets large. Figure|5]shows the transition of decoding 
error rate of 2D-SC {di = 3,dr = 6, L = 101, uj,Z) codes 
with w = 2 and Z as square segment of size 15. The channel is 
BEC(e = 0.48) with 20 burst section erasures injected. These 
burst section erasures are recovered by 2D-SC codes. It is 
observed that 2D-SC codes can recover more burst section 
erasures as L gets large. 

VI. Conclusion 

We propose MD-SC codes. We observed that 2D-SC codes 
are more robust to burst section erasures than ID-SC codes. 



0.49 
0.48 
0.47 
0.46 
0.45 
0.44 
0.43 
0.42 
0.41 
0.40 



Uncoupled MAP threshold 
2D Coupled BP Threshold 



10 



20 30 



40 50 

z 



60 70 



80 



90 100 



Fig. 3. The BP threshold of 2D-SC (d; =3,dr 
shortening domain Z as hypercube of size z. 
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Fig. 4. Transition of decoding error rate of 2D-SC (d; = 3, rfr = 6, L = 101, oj, Z) codes with w = 2 and Z as line segment of width 1. The channel is 
BEC(e = 0.48) with 20 burst section erasures injected. 




Fig. 5. Transition of decoding error rate of 2D-SC (d; = 3, dr = 6, L = 101, u), Z) codes with to = 2 and Z as square segment of size 15. The channel 
is BEC(e = 0.48) with 20 burst section erasures injected. 



